arXivrl 502.05045v3 [astro-ph.CO] 22 Mar 2016 


Mon. Not. R. Astron. Soc. 000, 1-?? (2014) 


Printed 24 March 2016 (MN IM’eN style file v2.2) 


Group-galaxy correlations in redshift space as a probe of 
the growth of structure 

F. G. Mohammad*^S. de la Torre^, D. BianchP’^, L. Guzzo^, J. A. Peacock^ 

^INAF-Osservatorio astronomico di Brera, via Emilio Bianchi 46, 1-23807 Merate (LC), Italy 
Universitd deWInsubria, Dipartimento di Scienza e Alta Tecnologia, Via Valleggio 11, 122100 Como, Italy 
^Aix-Marseille Universite, CNRS, LAM (Lahoratoire d'Astrophysique de Marseille) UMR 7326, 13388, Marseille, Franee 
"^ICG - Institute of Cosmology and Gravitation, University of Portsmouth, Dennis Seiama Building, Burnaby Road, 

Portsmouth, PCI 3FX, United Kingdom 

^Institute for Astronomy, University of Edinburgh, Royal Observatory, Blaekford Hill, Edinburgh EH9 3HJ, United Kingdom 


24 March 2016 


ABSTRACT 

We investigate the use of the cross-correlation between galaxies and galaxy groups to 
measure redshift-space distortions (RSD) and thus probe the growth rate of cosmo¬ 
logical structure. This is compared to the classical approach based on using galaxy 
auto-correlation. We make use of realistic simulated galaxy catalogues that have been 
constructed by populating simulated dark matter haloes with galaxies through halo 
occupation prescriptions. We adapt the classical RSD Dispersion model to the case 
of the group-galaxy cross-correlation function and estimate the RSD parameter (3 by 
fitting both the full anisotropic correlation function ^^(rp,7r) and its multipole mo¬ 
ments. In addition, we define a modified version of the latter statistics by truncating 
the multipole moments to exclude strongly non-linear distortions at small transverse 
scales. We fit these three observable quantities in our set of simulated galaxy cata¬ 
logues and estimate statistical and systematic errors on /3 for the case of galaxy-galaxy, 
group-group, and group-galaxy correlation functions. When ignoring off-diagonal el¬ 
ements of the covariance matrix in the fitting, the truncated multipole moments of 
the group-galaxy cross-correlation function provide the most accurate estimate, with 
systematic errors below 3% when fitting transverse scales larger than 10/i“^Mpc. In¬ 
cluding the full data covariance enlarges statistical errors but keep unchanged the level 
of systematic error. Although statistical errors are generally larger for groups, the use 
of group-galaxy cross-correlation can potentially allow the reduction of systematics 
while using simple linear or Dispersion models. 
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1 INTRODUCTION 

Cosmological observations over the past 15 years have es¬ 
tablished that the Universe is undergoing a late-time phase 
of accelerated expansion (e.g. Riess et al. 1998; Perlmutter 
et al. 1999). This is most naturally explained with the pres¬ 
ence of a ‘dark energy’, a fluid with negative pressure hll- 
ing the entire Universe, which is currently indistinguishable 
from Einstein’s cosmological constant. Alternatively, how¬ 
ever, one may consider reproducing observations by modify¬ 
ing the very nature of the gravitational equations of General 
Relativity (e.g. Carroll et al. 2004; Joyce et al. 2014). Geo¬ 
metrical probes such as the Cosmic Microwave Background 
(CMB), Baryon Acoustic Oscillations (BAO) and type la 
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Supernovae (SNla) constrain the expansion history H[z), 
which however can be equally well fitted by both scenarios. 
This degeneracy between dark energy and a modihcation 
of standard gravity can only be lifted by dynamical probes 
looking at the growth of structure inside the Universe, which 
is directly linked to the underlying theory of gravity (Joyce 
et al. 2014). 

The growth of cosmological structure induces galaxy pe¬ 
culiar velocities, i.e. coherent flows of galaxies towards mat¬ 
ter overdensities. When redshifts are used to map galaxy 
positions, we are sensitive to such peculiar velocities, whose 
line-of-sight component combines with the cosmological red- 
shift. As a result, the reconstructed spatial distribution of 
objects is distorted in the radial direction, what are referred 
to as redshift-space distortions (RSD). RSD can be quan¬ 
tified statistically in galaxy redshift surveys by modelling 
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the corresponding anisotropy that can be measured in the 
two-point correlation function (2PCF) ^(r) or, correspond¬ 
ingly, in the power spectrum P(k). In the linear regime, 
this effect is described in Fourier space by the Kaiser model 
(Kaiser 1987), or its equivalent in configuration space devel¬ 
oped by Hamilton (1992). Such modelling is however com¬ 
plicated by the non-linear evolution of matter overdensities 
and velocities (see e.g. de la Torre & Guzzo 2012, and refer¬ 
ences therein). 

On large scales we can still observe the linear, coherent 
flows tracing the growth of cosmological structure, which 
enhance the amplitude of two-point correlations. But small- 
scale motions are dominated by high-velocity galaxies in viri- 
alized structures, resulting from the dynamical evolution of 
the highest peaks in the density field. The resulting stretch¬ 
ing effect in galaxy maps is commonly referred to as the 
Fingers of God effect (FoG), due to the elongated shapes 
that groups and clusters acquire when observed in redshift 
space. These two regimes of RSD produce characteristic fea¬ 
tures in the two-point statistics of the galaxy distribution, 
which can be studied by computing the correlation function 
as a function of two variables, Vp and tt, perpendicular and 
parallel to the line of sight, respectively. The two effects, 
linear and non-linear, introduce respectively a large-scale 
squashing and a small-scale elongation along the line-of-sight 
direction tt. In linear theory, the amplitude of the squashing 
effect is directly proportional to the logarithmic growth rate 
of density fluctuations f{z) (see Eq. (8)) (Kaiser 1987). In 
practice, our discrete tracers will generally be biased with 
respect to the overall matter distribution. Under the hypoth¬ 
esis of a linear bias biin{z), the large-scale squashing effect 
will depend on the parameter /3(z) = f{z)/b\[n{z). 

An empirical correction to the linear model was intro¬ 
duced as to account for the contribution of nonlinear distor¬ 
tions, what is normally referred to as the Dispersion model 
(Fisher et al. 1994; Peacock & Dodds 1994). In this model, 
the linear expression is convolved (in configuration space), 
or multiplied (in Fourier space) with a given distribution 
function for the pairwise velocities (p{v\\) along the line of 
sight. The dispersion model has been widely and success¬ 
fully used in the past to estimate /3(z) from measurements 
of ^(rp,7r) (Peacock et al. 2001; Hawkins et al. 2003; Guzzo 
et al. 2008; Ross et al. 2007; Gabre & Gaztahaga 2009a,b; 
Gontreras et al. 2013) or the power spectrum (Percival et 
al. 2004; Tegmark et al. 2004; Blake et al. 2011). It has be¬ 
come clear that such empirical correction provides estimates 
that can be biased by up to 10%, typically depending on the 
bias of the tracer being employed (Okumura & Jing 2011; 
Bianchi et al. 2012). The precision already achieved by cur¬ 
rent largest surveys requires to improve on these limitations 
to reach systematic errors of the order of the percent (e.g. 
BOSS, Reid et al. 2014; Rentier et al. 2014; Samushia et al. 
2014; Sanchez et al. 2014; Alam et al. 2015). Such limita¬ 
tions are intrinsic in the empirical nature of the dispersion 
model (Scoccimarro 2004; de la Torre & Guzzo 2012) and 
significant effort has been invested over the past five years to 
improve this theoretical description (e.g. Scoccimarro 2004; 
Taruya et al. 2010; Reid & White 2011; de la Torre & Guzzo 
2012, and references therein, Reid et al. 2014). Alternatively, 
these theoretical limitations may be evaded if we are able 
to (a) identify galaxy tracers that are less affected by non- 
linearities, and/or (b) build novel statistics of galaxy cluster¬ 


ing for which model fits are less sensitive to the same effect, 
while keeping the theoretical model as simple as possible. 

In this paper we investigate two options along these two 
avenues. Specifically, we consider the case in which one has 
access to a robust galaxy group catalogue within the vol¬ 
ume of a corresponding galaxy redshift survey; this will be 
possible with next-generation galaxy redshift surveys such as 
Euclid (Laureijs et al. 2011), also in combination with obser¬ 
vations in other bands, as in particular the X-ray data from 
e-Rosita (Merloni et al. 2012). We therefore first consider the 
merits of using the galaxy group-galaxy cross-correlation to 
extract the distortion parameter [3. We might expect the 
centre of mass of a galaxy group to have little or no ran¬ 
dom velocity, since the internal orbital velocities produce 
Eingers-of-God only in the overall galaxy distribution. In 
this view, group-group correlations should be ideally placed 
(as also suggested by the dependence of the systematic er¬ 
ror on the halo mass shown by Bianchi et al. 2012). More¬ 
over, the group-galaxy cross-correlation provides a higher 
statistical power than using group auto-correlation, given 
the higher number density of the galaxy catalogue. Secondly, 
we investigate a new way of compressing the cosmological 
information present in the anisotropic two-point correlation 
function in redshift space, by building modified (truncated) 
multipole moments of ^(rp,7r) that partially mitigate small- 
scale non-linearities. These functions avoid explicitly includ¬ 
ing the contributions from small transverse scales, domi¬ 
nated by non-linear distortions^. In this analysis, we there¬ 
fore compare the standard approach of using galaxy or group 
auto-correlation and the novel approach of using the galaxy- 
group cross-correlation, as well as study how the different es¬ 
timators of the two-point statistics, namely the anisotropic 
2PGE ^^(rp, tt), its standard multipoles ^^’^^^(s) and its trun¬ 
cated multipoles (s), behave in recovering the distortion 
parameter [3 in simulations. 

The paper is organized as follows. In Section 2, we de¬ 
rive the linear Kaiser/Hamilton and Dispersion models for 
the cross-correlation in the case of different observables: the 
anisotropic 2PGE ^^(rp,7r), its multipole moments 
and the related truncated multipole moments ^^’^^^(s). In 
section 3, we describe the data, two-point correlation func¬ 
tion estimators, as well as ingredients needed to build the¬ 
oretical models. In Section 4 we present the comparative 
analysis of the different approaches and methods. Einally, 
in Section 5 we discuss our results and conclude. 


2 REDSHIFT-SPACE DISTORTIONS: 
MODELLING 

The apparent positions of objects are modified if we use 
redshifts to infer cosmological distances. The line-of-sight 
component of the peculiar velocity distorts positions in the 
following way: 

s = r + u(r)-e||, (1) 

where s and r are objects positions, respectively, in redshift- 
and real space; e|| is the unit vector along the line of sight; 

^ In a parallel work, Reid et al. 2014 independently define a sim¬ 
ilar statistics as to mitigate fibre collisions effects in the analysis 
of RSD in the BOSS survey. 
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u is the scaled velocity field defined as u(r) = v(r)/aif(a) 
with v(r) being the peculiar velocity field; a is the scale 
factor; and H(a) is the expansion rate of the Universe. The 
s coordinates constitute the so-called redshift space and the 
distortions produced in the matter distribution are usually 
referred to as Redshift-Space Distortions (RSD). 

2.1 Linear Model for the Cross-Correlation 

We first derive the linear model for the two-point cross¬ 
correlation function. Following the derivation of Kaiser 
(1987) for the auto-correlation, we start with assuming mass 
conservation in real- and redshift space in terms of overden¬ 
sities (5(r) = p(r)/(p) — 1: 

[1 + (5m(s)] = [1 + (5m(r)] (2) 

where (5m (r) and (5m (s) are, respectively, the matter density 
contrast in real- and redshift space. The volume element 
in redshift space df's is related to that in real space 
through d^s — |J|d^r where \J\ is the Jacobian of the re¬ 
lated coordinates transformation. In the particular case of 
the transformation in equation (1), the Jacobian is given by 

|J| = l + a||M||(r), (3) 

where d\\ = d/dr\\ and r^||(r) = u(r) • e|| is the line-of-sight 
component of the scaled velocity field u(r). Equation (2) can 
thus be written as 

5m(s) = [(5,„(r)-a||M||(r)] [1+ a||M||(r)]“\ (4) 

In the linear regime approximation (i.e. small real-space 
density field |(5(r)| ^ 1, small velocity gradients |^||'?^|| | ^ 1 
and an irrotational velocity held Vxu(r) = 0), equation (4) 
becomes 

<5^(8)= (r) - a, I W| I (r)] [l - S,, w,, (r)] . (5) 

Neglecting higher order terms in equation (5) gives 

5m(s) = 5m(r) - a||M||(r). (6) 

In the plane parallel approximation, when the scales of per¬ 
turbations are assumed to be much smaller then their dis¬ 
tances from us (i.e. u\\/r\\{r) ^ I), equation (6) can be 
written in Fourier space as 

SL{k,iJ,k) = [l + f{z)ul] 5m.{k). (7) 

with Ilk = k\\/k. In equation (7) / is the growth rate of 
structure dehned as: 

dln^ (8) 

Equation (7) is valid for the overall matter density held. 
However, when objects of a given type i are used as tracers 
of the overall matter density held, their distribution is in 
general biased with respect to the underlying matter. At 
hrst order, the overdensity 5i of a given population i can be 
related to that of the overall matter 5m through a linear bias 
factor bi 

(54k) = 5, (5m (k). (9) 

In this work we assume the linear bias bi to be a scale- 
independent parameter which is valid with good approxi¬ 
mation in regimes that are not affected by small-scale non 


linearities and the ones due to the Baryon Acoustic Oscilla¬ 
tions (BAO). 

We also assume the peculiar velocity held of z-type ob¬ 
jects to be unbiased with respect to that of the overall mat¬ 
ter. This assumption holds on scales much larger than the 
typical size of virialized structures. Indeed, inside such struc¬ 
tures dynamical processes such as dynamical friction and 
tidal disruption may alter the velocity held of objects i with 
respect to that of the total matter introducing a further 
velocity bias (Munari et al. 2013; Elia et al. 2012). Thus, 
equation (7) becomes 

5l{k,Uk) = Si{k) + f{z)n\5m{k). (10) 

Dehning the distortion parameter 44^), related to the 
population 4 as 

= si' 

and combining it with equation (10) and the linear bias re¬ 
lation (9) leads to 

5l{k,iik) = [l l3i{z)nV\ 5i{k). (12) 

Given two population of objects, e.g. individual galaxies 
‘gaV and their groups ‘gr\ their cross power spectrum is 
dehned as 

(5gal(kl)<5g,(k2)) = (27r)"fo(k2 - kl)Pcr(k). (13) 

The linear cross power spectrum in redshift space A1in(k) 
can be related to that in real space Pcr(k) by combining 
equation (13) with (12) 

f’cr,lin(fc, Uk) = [1 + /^galMfc] [l + hrul] -Pcr(fc). (14) 

It is useful to write Pcr,iin(^? Mfc) as a sum of spherical har¬ 
monics 

P!rMk, Hk) = E P:ixlik)LdUk), (15) 

here, Li(jik) are the Legendre polynomials and ■IxSk) are 
the multipole moments of the linear cross power spectrum 
■Pcr,iin(k) given by 

P:ifnik) = ^4^ 1^' PdMk,Uk)Le{pk)duk. (16) 

The equivalent expression for the two-point cross¬ 

correlation function ^4,iin(^^) is provided by Hamilton (1992) 

Ccr.lin(rp,7r) = (17) 

£ 

In equation (17), and tt are respectively the components 
of the pair separation s transverse and parallel to the line 
of sight, pi is the cosine of the angle between the pair sepa¬ 
ration s and the line of sight and multipole 

moments of ^4,iin(^p?' tt): 

e!i?n(«) = *4 '^p:;!Dk)Mks), (18) 

where ji{ks) are the spherical Bessel functions. The only 
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non-null multipole moments are 








1+ gd^sal + /5gr) + g/5gal/5gr 


'2 4 

~ (^gal + /^gr) + — /^gal/^gr 


^/3gal/3gr 


^cr(r) (19a) 
[6r(r)-^cr(r)] (19b) 

(19c) 


^cr(^)+ 2^cr(^) 2^cr(^) 


where ^cr(r’) and ^cr(r’) are the integrals of the real-space 
angle-averaged two-point cross-correlation function ^cr{r) 
(Cole et ah 1994): 

^cr(r) = 4" / ^cr{r')r'‘^dr' (20a) 

lcr(r) = [ Gddy^dr'■ (20b) 

^ Jo 

This model was adopted by Mountrichas et ah (2009) to 
measure the bias of the QSOs in 2SLAQ, 2QZ and SDSS, 
by cross-correlating them with a population of Luminous 
Red Galaxies (LRG) in 2SLAQ and AAOmega. 

Both /3gai and ^dgr in fact describe the same growth rate 
/ (z); therefore using the linear bias relation (9), the 2PCF 
of the population z, in real space, can be written as 

Ci(r) = 6iCm(r). (21) 


Combining equation (21), written separately for galaxies ^gai 
and groups ^gr, with (11) gives 


/3gr(2) = &12(2:)/3gal(«). (22) 


In equation (22), bi 2 {z) is the ‘relative bias’ between galaxies 
and groups, dehned as 


bi2{z) 



&gal(g) 

6gr(z) 


(23) 


This quantity can be estimated directly from the data, 
once the real-space correlation functions of the two popula¬ 
tions are measured through projection of <f^(r’p,7r). Dehned 
in this way, bi 2 will be smaller than unity, given the larger 
bias of groups with respect to galaxies. We prefer this def¬ 
inition as it allows us to obtain a more compact notation 
in the following equations. Using equation (22), therefore, 
equations (19a),(19b), (19c) become 


C4°l(s) = 

gJZg) = 

&s) = 


1 + g/^ga^l + 612 ) + g/5|al&12 


2/5gal(l + &I 2 ) +;^^gal&12 


2g/5gal&12 


Gr{r) + ^Gr{r)-'^Gr{r) 


fcr(r) (24a) 
[ecr(r)-^„(r)] (24b) 
• (24c) 


The linear Kaiser/Hamilton model for the two-point 
auto-correlation function is recovered just by taking bi 2 {z) = 
1 and replacing the real-space two-point cross-correlation 
function Ccr{r) and its integrals ^cr{r), ^cr{r) with their coun¬ 
terparts in the auto-correlation case. 


2.2 The Dispersion Model 

To derive the linear model in Section 2.1, strong assumptions 
have been made, limiting its validity to the very large scales. 
We adopt the mostly used empirical ‘Dispersion model’ 


(Peacock & Dodds 1994) to model the shape of the 2PGF on 
small and intermediate trans-linear scales. This model con¬ 
volves the linear model presented in Section 2.1 with a pe¬ 
culiar Pairwise Velocity Distribution (PVD) function along 
the line of sight 


G{rp,n) = 


-f 


Cll. 


rp,TT - 


^’11 

aH{a) 


ip{vGdv\\. 


(25) 


Different functional forms have been proposed in the litera¬ 
ture for the PVD function. We adopt an exponential distri¬ 
bution function for which is found to be in good agree¬ 

ment both with data from N-body simulations and with ob¬ 
servations from large galaxy redshift surveys (Zurek et al. 
1994): 

1 ^hlll 

=-e -12 . (26) 

V2cri2 

Here, cri 2 is usually referred to as the peculiar pairwise ve¬ 
locity dispersion along the line of sight. In this paper we 
assume <712 to be a scale-independent free parameter. More 
complicated and accurate models for the peculiar pairwise 
velocity distribution have been calibrated on simulations to 
take into account the scale dependence of the velocity dis¬ 
tribution (e.g. Zu V Weinberg 2013; Bianchi et al. 2014). 


2.3 Truncated Multipole Moments 

The standard multipole moments of the 2PGF ^^(s) are ob¬ 
tained by projecting it onto the Legendre polynomials: 

G{s,FMn)dn, (27) 

with = Tp + TT^ and /i = tt/s. 

Since strong non-linear distortions affect small trans¬ 
verse (rp) scales, we propose an alternative way of using the 
multipole moments of the 2PGF by removing such scales. 
This means that, at a given scale s, we consider only the 
clustering signal at rp > fp. The new ‘multipole moments’ 
are then given by 

^ J G{s,FLe{ydfj., (28) 

where 

It is important to stress here that although ^^’^^^(s) are 
not, mathematically speaking, the multipole moments of the 
2PGF, in the following part, for the sake of simplicity, we 
will refer to them as the ‘truneated multipole moments’. 

The gain in using the truncated multipole moments 
is given by the fact that they reconcile the classical ap¬ 
proach of using the anisotropic 2PGF ^(rp,7r) with the 
one of using its multipole moments. Indeed, the anisotropic 
2PGF ^(rp,7r) allows us to systematically exclude from the 
ht small transverse rp scales, affected by strong non-linear 
distortions, which are difficult to model analytically. How¬ 
ever, this approach would involve a huge covariance ma¬ 
trix: it is thus in practice computationally infeasible. The 
size of the covariance matrix in the case of multipole mo¬ 
ments is much smaller, but nonlinearities on small transverse 
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Figure 1. Real-space two-point correlation function measurements from the MultiDark Runl data at 2 ; = 0.1, averaged over 27 equal 
sub-samples. Blue line shows the (HOD, L > IL*) galaxy auto-correlation, red line represents the auto-correlation of (M > Mq 

dark matter haloes) groups while the group-galaxy cross-correlation is plotted with green line. The CAMB model prediction for the 
correlation function of the underlying matter density field is also plotted (black line). The red dashed line is the power-law extrapolation 
of the group auto-correlation function on scales below 3 h“^Mpc using parameters in Table 1. The relative statical errors on the averaged 
correlation functions are shown as the transparent filled contours. 


scales make an undesired contribution when projecting the 
redshift-space 2PCF ^^(s, /i) onto the Legendre polynomials. 
The truncated multipoles allow us to remove these strong 
non-linear distortions, while retaining a relatively small co- 
variance matrix and yielding a numerically tractable prob¬ 
lem. 


3 METHODOLOGY 

To test the model presented in Section 2, we make use 
of simulated catalogues of galaxies and groups based on 
the MultiDark Runl (MDRl) dark matter N-body simu¬ 
lation (Prada et al. 2012). MDRl is a high mass resolu¬ 
tion simulation within a cube of side 1000/i“^Mpc, with 
Np = 2048^ dark matter particles and the mass of each par¬ 
ticle being Mp = 8.721 x 10^Mq. It assumes a ACDM 
cosmology using cosmological parameters from WMAP5 
and WMAP7 data releases (Qrn, ^6 ,'/T's,/ i, c^s) = 
(0.27,0.73,0.0469,0.95,0.7,0.82). Dark matter haloes in 
MDRl are identified using the Friends-of-Friends (FoF) al¬ 
gorithm. 

We use the simulated galaxy catalogue used in de 
la Torre & Guzzo 2012, which was constructed by pop¬ 


ulating MDRl dark matter haloes and specifying the 
Halo Occupation Distribution (HOD). It contains 2,945,687 
galaxies with a B-band luminosity of L > L*. We 

used the HOD parametrization of Zheng et al. (2005) 
with parameters (log Mmin, cr\ogm Aog Mq, log Mi, a) = 
(12.18,0.21,12.18,13.31,1.08) with masses given in units of 
MqIH. The details of the HOD model that has been used 
and the construction of the sample are provided in Appendix 
B of de la Torre & Guzzo 2012. We only remark here that ha¬ 
los are assumed isotropic and spherical, with satellite galaxy 
velocities drawn from Gaussian distribution functions along 
each Gartesian direction, with dispersion computed follow¬ 
ing van den Bosch et al. (2004). 

For the group catalogue, we use all dark matter haloes 
with mass M > 10^^Mq giving a total of 350,518 groups. 
This specific mass threshold is chosen as to match the ob¬ 
served number density of groups in the 2dFGRS Percolation- 
Inferred Galaxy Group (2PIGG) catalogue drawn from the 
2dFGRS (Eke et al. 2004). 

3.1 Estimator 

We measure the 2PGF using the minimum variance estima¬ 
tor proposed by Landy & Szalay (1993) adapted to the case 















6 F. G. Mohammad et. al. 


of cross-correlation function: 

, DiD 2 {s) - DiR{s) - D 2 R{s) + RR{s) , 

= - RRG) -’ 

where D 1 D 2 are the data — data pairs between the two dif¬ 
ferent catalogues, DiR are the datai — random pairs and 


RR are the random — random pairs. All pair counts are 
normalized to the related total numbers of distinct pairs. 

Multipole moments of the measured two-point correla¬ 
tion function ^(s) are obtained by simply projecting it onto 
the Legendre polynomials (equation 31): 




2£+l 


/:[ 


DiD 2 (s, jj) — DiR{s, jj) — D 2 R{s, jj) + RR{s, ji) 


RR{s, ji) 


Li{ia)dia 


r 


(31) 


To measure the truncated multipole moments, contain¬ 
ing the clustering signal only on transverse scales > fp, 
the integration in equation (31) is truncated at [—/x, +/x] 
rather than [—1, +1] with fl given by equation (29). We make 
use of a random catalogue with a number of random points 
10 times larger than the galaxy sample to reduce the shot 
noise in our measurements. 

We use logarithmic s bins of size Asiog = 0.02 that 
covers the range [0.1,200] h~^Mpc. fi is divided into 200 
linear bins between [0,1]. The measurements of the real- 
space correlation functions are shown in Figure 1, while the 
measured standard and truncated multipoles for the case of 
fp = 5/i“^Mpc are shown in Figure 2. For the anisotropic 
2PCF ^(rp,7r), both transverse rp and parallel to the line- 
of-sight TT components of pair separation are linearly binned 
with bins of size 1 /i“^Mpc in the interval [0,100] /i“^Mpc. 
The measured anisotropic 2PCF are presented in Figure 3. 

loglo(»'i+l) = loglo(»'i) + Slog (32) 

Measurements in linear bins are sampled at the mid point of 
each bin while in the case of logarithmic binning a logarith¬ 
mic average (eq. 33) of the two edges of each bin is taken as 
the reference point. 

In the last part of the analysis, we fix the number of bins in 
s to Nb = 12 inside the range [smin, 80]/i“^Mpc in order to 
fit multipole moments with the full covariance matrix (see 
Section 3.2). 

3.2 Error Covariance Matrix 

In principle, the determination of the covariance matrix 
requires many independent realizations of our dataset, of 
which we inevitably possess only a single example. As an 
approximation to the ideal case, essentially two classes of 
methods have been proposed: (a) use simulated data, to pro¬ 
duce Nr independent mock samples with properties as close 
as possible to the observed data; (b) use internal methods, 
in which Nr multiple realizations are constructed from the 
overall data, by resampling the observations in some appro¬ 
priate way. In particular, four such methods have been used 
in the past literature: (1) the classical bootstrap method, in 
which each realization is constructed by resampling with 
replacement of single objects in the data set. In the other 
three methods, the sample volume is split into a number 
of sub-volumes, which are then combined following different 


schemes: (2) the bloek-wise bootstrap method builds realiza¬ 
tions by combining all sub-volumes, but assigning a random 
weight to each of them, to test the sensitivity of the mea¬ 
sured statistics with respect to specific parts of the sam¬ 
ple; (3) the jaek-knife method builds realizations simply ob¬ 
tained by omitting one of the sub-volumes; finally (4) the 
sub-sample method treats each sub-volume as a realization 
of the available dataset, yielding a single estimate of the set 
of the statistics under evaluation. 

Of these internal methods, the sub-sample method is 
distinct, as it explicitly considers datasets whose volumes 
are lower than the original parent sample; conversely, the 
bootstrap and jack-knife methods attempt to estimate a co- 
variance matrix that is appropriate for the whole sample. 
However, if we are willing to work with smaller volumes, the 
sub-sample method is closer to the ideal case of many re¬ 
alizations. The sub-samples are not truly independent, and 
so the results will not be correct on scales approaching that 
of the sub-volume, but otherwise the sub-sample method 
should yield directly a reliable covariance matrix for samples 
having the size of a single sub-volume. If we divide the ini¬ 
tial dataset into Ns sub-volumes, we can derive the desired 
statistic (anisotropic 2PCF or its multipoles, in this case) 
for each sub-volume, so the natural thing to do is to average 
these sub-estimates (as opposed to estimating the statistic 
over the whole sample). In the small-scale limit where the 
sub-samples can be treated as independent, the covariance 
matrix for this mean statistic would then be just reduced 
by a factor Ns, since covariances of independent samples 
add linearly. The overall covariance matrix for the 2PCF, 
measured as an average over sub-samples, is therefore 

^ 

Here, yf is the measurement of the statistic y in bin i made 
from the k-th realization while {yi) is the mean value of y in 
the same bin over the Ns realizations. We emphasise again 
that the quantity {yi) is to be used as the result for analysis, 
and it is not identical to the same statistic evaluated once 
over the whole volume. In practice, we divided the dataset 
into Ns — 3^ sub-cubes, so that the side of a single sub¬ 
cube is 333/i“^Mpc. This is large enough that the scales of 
interest for RSD (tens of Mpc) can be treated as independent 
in each sub-volume. 

Clustering measurements show strong bin-to-bin cor¬ 
relations that needs to be taken into account in the fit¬ 
ting procedure using the covariance matrix. Because the 
anisotropic 2PCF has a large number of separation bins. 
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(a) Standard multipoles 



(b) Truncated multipoles 

Figure 2. Standard (top panel) and truncated (bottom panel) monopole i = 0 and quadrupole ^ = 2 of the 2PCF in redshift space 
averaged among Ng = 27 equal sub-volumes drawn from the MDRl simulation at 2 : = 0.1. The blue dashed lines show the auto correlation 
of (HOD, L > IL* ) galaxies, red dash-dotted lines represent the group (M > 10^^h~^M q dark matter haloes) auto-correlation and 
their cross correlation is plotted with green continuous lines. The shaded regions show the related rms among Ng measurements scaled 
by 1/Ng. 






















8 F. G. Mohammad et. al. 




Tp [h ^Mpc] 


(c) Cross-correlation 

Figure 3. Anisotropic two-point auto-correlation functions C^(rp,7r) for galaxies (top left panel), groups (top right panel) and their 
cross-correlation (bottom panel) averaged among Ng = 27 equal sub-volumes. 


the measurement of a proper covariance matrix requires a 
huge number of independent realizations which is not avail¬ 
able here. Therefore, in order to fairly compare [3 parame¬ 
ters obtained from the anisotropic 2PCF and its multipole 
moments, we first restricted the analysis to using diagonal 
errors. We then used the full covariance matrix but only 
to compare results obtained from multipole moments. The 
multipoles have a smaller number of bins which makes eas¬ 
ier the estimate of the full covariance matrix given a limited 
number of realisations. In general, the number of bins in the 
multipole moments of the 2PCF that can be used should 
be smaller than the number of independent realisations. In¬ 
deed, a higher number of bins would inevitably yield to a 
singular covariance matrix (e.g. Hartlap et al. 2007) which 
cannot be inverted. To estimate the covariance matrix we 
thus fixed the number of bins to be iV;, = 12 on each multi¬ 
pole, independently of the fitting range we consider. In the 


case of the truncated multipoles, however, the measurement 
in the first bin is always equal to zero. To overcome this issue 
we exclude the first bin from the fitting procedure and esti¬ 
mate the covariance matrix on the remaining bins resulting 
in iVfo = 11. 

Given the small number of independent realizations 
{Ng = 27) that we possess and that is of same order as 
the number of bins in our measurements, the raw esti¬ 
mate of the covariance matrix is not well constrained and 
rather noisy. Therefore, we adopt the method of Mandel- 
baum et al. 2013 to reduce its noise, by smoothing the 
off-diagonal elements of the associated correlation matrix 
Rij = Cij with a box-car algorithm. For this, 
we use an optimal box size of 3 x 3 that avoids altering 
the correlation matrix structure. It is important to stress 
that in this procedure, the smoothing has been applied sepa¬ 
rately for each sub-quadrant and the diagonal elements have 
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(a) Covariance matrix 


(b) Correlation matrix 


Figure 4. Covariance (left panel) and correlation (right panel) matrices of the truncated multipoles of the two-point group-galaxy 
cross-correlation between [5, 80]/i“^Mpc. In each panel the matrix above the white continuous line shows the direct measurement from 
27 equal sub-volumes while the lower part show the matrix after it has been smoothed using a box-car algorithm on 3 x 3 sub-matrices 
(right panel). 


been left unchanged. As an illustration, we show in Figure 
4 both the original and smoothed covariance and correla¬ 
tion matrices, in the case of the truncated multipoles of the 
group-galaxy two-point cross-correlation function in range 
s = [5, 80]/i“^Mpc. One can see in this figure how the global 
structure of the original covariance matrix is preserved in 
by smoothing process and the noise is reduced for the very 
off-diagonal elements. 


-21nr(^,(7) = x^(/3,o-) = [y*''* 


here, C is the likelihood function, and are 

respectively the measurement and the model prediction for 
the htted quantity in bin i and C~-^ is the inverse covariance 
matrix. 

In particular, since our model depends only on two 
parameters (/3,cr), we explore the parameter space using 
a grid on /3 G [0,0.60] with bins of size A/3 = 10“^ and 
cfi 2 G [0,10] /i“^Mpc with bins of Acri 2 = 0.05/i“^Mpc. 


^ We will express the dispersion parameter in units of length 
[h“^Mpc]. 


3.3 Fitting Method 

The model we presented in Section 2 depends mainly on 
two free parameters: the distortion parameter and the 
dispersion parameter cri 2 ^. We ht the measured two-point 
statistics in redshift space with their models (Section 2) by 
minimizing the quantity in equation (35): 


< 7 )] • C-f ■ [j/f * - yr\0, < 7 )] (35) 


3.4 Model Construction 

In this section we present the measurements, from the sim¬ 
ulated catalogues, of the main ingredients required to con¬ 
struct the model presented in Section 2. 


3.4-1 Real-Space Correlation 

One such ingredient is the angle-averaged 2PCF in real space 
^(r). Given our aim here, which is to test the relative perfor¬ 
mances of different estimators of the growth rate, we use the 
result directly available from the simulation itself. For real 
data this quantity is also available in principle, either via 
the projection of ^^(r’p,7r), or by htting an assumed model. 
Specihcally, as noted above, we use the real-space correlation 
functions averaged among 27 sub-sample realizations as the 
input for the model that we denote simply with ^(r) in the 
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Figure 5. Power-law fits of the real-space two-point correlation 
functions, averaged among 27 equal sub-samples, for galaxy auto¬ 
correlation (blue dashed line and squares) and the group-galaxy 
cross-correlation (green continuous lines and circles). Points with 
error bars represent the measurements while the lines result from 
the best fit of data using equation (36) (Table 1). 
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Figure 6. Same as in Figure 5 but for the case of the group 
auto-correlation function. 


following. In Figure 1 we show the measurements of such 
two-point correlation functions for galaxy auto-correlation 
(blue line), group auto-correlation (red line) and for the 
group-galaxy cross-correlation (green line). The correlation 
function for the underlying overall matter density field as 
predicted by the CAMB model (Lewis et al. 2000) is also 
presented (black line). We plot the quantity r‘^^{r) rather 
than simply ^(r), to enhance the differences between differ¬ 
ent correlation functions. The transparent filled areas rep¬ 
resent the statistical errors on the averaged real-space two- 
point correlation functions. 


3.4-2 Power-Law Extrapolation 


Although the linear model at a given scale r depends on 
the integrals of the real-space angle-averaged correlation 
function ^(r) between separation zero and r, our measure¬ 
ments of ^(r) are performed within [0.1,200] /i“^Mpc. We 
use a power-law form to extrapolate ^(r) on scales below 
0.1/i“^Mpc. 


^(r < O.l/i ^Mpc) 



(36) 


The power-law parameters (ro,7) are measured by fitting 
^(r) on very small scales, i.e. [0.1,0.25] h~^Mpc in the case 



ro [h ^Mpc] 

7 

Galaxy auto-correlation 

o q + 0.7 

, 7^ + 0.03 
' ^-0.03 

Group auto-correlation 

o rq + 0.02 

».py_o 02 

1 c;eq“l“0-006 

j..oo»_o 006 

Group-galaxy cross-correlation 

5 6 + °-3 

q q 1 -|-0 . 03 
^•^-L_o.03 


Table 1. Power-law ^(r) = (r/ro)~^ parameters estimated fit¬ 
ting the real-space two-point auto-correlation function of galax¬ 
ies and the group-galaxy cross-correlation function, averaged over 
27 equal sub-samples, between [0.1,0.25] h“^Mpc. The same pa¬ 
rameters in the case of the group auto-correlation function are 
estimated fitting it between [3,10] h“^Mpc. 



Figure 7. Measurements of the galaxy linear bias 5gal(r) and 
related l-cr errors are shown with blue dashed line and filled blue 
contours. The red continuous line with red filled contours rep¬ 
resent the measurement of the groups linear bias bgr{r) and its 
l-fj statistical errors. The black dashed and continuous lines with 
grey contours represent the best fit of data, respectively for the 
galaxy and groups linear bias factor, between [10,50] h“^Mpc 
with a constant model = const. 


of the galaxy auto-correlation and the group-galaxy cross¬ 
correlation functions. For the group auto-correlation, given 
the low clustering due to the low number pairs of group¬ 
sized dark matter haloes on scales below ~ l-2h~^ Mpc, 
we extrapolate ^gvir) on scales below 3/i~^Mpc and the 
power-law parameters are obtained fitting the measured ^gr 
in [3,10] h~^Mpc. The fitting results are listed in Table 1 
and shown in Figure 5 for the cross-correlation (green line) 
and the galaxy auto-correlation (blue line) and in Figure 6 
for the group auto-correlation (red line). In Figures 5 and 
6, the error bars related to the statistical errors on the mea¬ 
sured correlation functions are smaller than the size of the 
plot symbols. 


3.4-3 ^(r) Integrals 

Here we compute the integrals of the real-space angle- 
averaged correlation functions ^(r) defined in equations 
(20a) and (20b). The integration is performed by splitting 
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Figure 8. Systematic errors on the galaxy distortion parameter estimated by fitting the anisotropic (HOD, L > IL*) galaxy auto¬ 
correlation function (blue lines) and the group-galaxy cross-correlation function (green lines) and on the group distortion parameter /3gr 
from the anisotropic {M > Mq dark matter haloes) group auto-correlation function (red lines), using only the diagonal elements 

of the covariance matrix, are plotted. The continuous lines with filled points represent results using the Dispersion model. The dashed 
lines with empty points correspond to fits using the linear Kaiser/Hamilton model only. The error bars correspond to the scatter among 
the 27 sub-samples used in each analysis. The right panel is just a zoom of the left panel. The shaded regions represent the 5% and 10% 
levels in the left panel while in the right panel are shifted to 3% and 5% level. 


it into two parts: 




= )l 1 


f r'‘^^{r)dr' 

J 

[ (37a) 

^ \ 


f r'^^(r)dr 

J 

[ (37b) 


Here rpi corresponds to the upper limit of the power-law 
extrapolation, which is rpi = 0.1/i~^Mpc for the galaxy auto¬ 
correlation and the group-galaxy cross-correlation and rpi = 
3/i“^Mpc for the group auto-correlation. In (37a) and (37b) 
the first terms come from the power-laws on scales between 
[0,rpi] /i“^Mpc, while the second terms are the integrals of 
the effective measurements of the real-space 2PCFs. 


3.5 Fiducial Model 


Among others, one of the advantages in using data from sim¬ 
ulations is that we know a priori the expected (i.e. fiducial) 
values of the quantities we want to measure. This makes 
them an ideal tool to test and compare the reliability of 
different methods and theoretical models. In this part we 
discuss how the fiducial values of the galaxy and group dis¬ 
tortion parameters, /3gai(^ = 0.1) and ^gr(« = 0.1) have been 
estimated from the simulated datasets. 

To perform this, we need first to estimate the fidu¬ 
cial value for the linear bias bi of our galaxies and groups 
(haloes). We do this by fitting the quantity 


bi{r) = 




1 1/2 


_ <^mass (l") _ 


(38) 


in equation (38) with a constant between [10,50] h~^Mpc. 
This choice for the fitting range is dictated by the need to 
exclude both small non-linear scales and BAO scales. Not 
having at our disposal the catalogue of dark matter particles 
for the MDRl, we have used CAMB to predict the matter 
correlation function ^mass(i") at z = 0.1. The ratio in equa¬ 
tion (38) is shown in Figure 7 as measured for the galaxies 
(blue dashed line and contours) and groups (red continuous 


line and contours). We obtain the following values: 

btt\iz = 0.1) = (1.331 ± 0.004) (39a) 

blf{z = 0.1) = (1.588 ± 0.003) (39b) 

The loss of power on scales > 30/i“^Mpc is not problematic 
since, given the statistical error bars, it is still consistent 
with the assumption of a constant bias factor. 

Once the two values of the linear bias hgai and bgr are 
measured, the estimation of the relative bias bi 2 is straight¬ 
forward using eq. (23) and results in: 

bi 2 iz = 0.1) = (0.838 zb 0.003) (40) 

Given the cosmology of the MDRl simulation, we ob¬ 
tain a fiducial value f^^{z — 0.1) = 0.5434, corresponding to 
the following values for the distortion parameters /3gai and 
/3gr: 

I3gii{z = 0.1) = (0.4083 ± 0.0006) (41a) 

,5“(« = 0.1) = (0.3422 ± 0.0006) (41b) 


4 RESULTS 

We now present the results obtained by fitting the measured 
two-point correlation functions with the models presented in 
Section 2. We shall first fit the full ^^(rp,7r), its multipoles 
^s,(£)(^) related truncated multipoles ^^’^^^(s) using 

only the diagonal elements of the data covariance matrices; 
we then perform a full covariance analysis using the trun¬ 
cated multipoles ^^’^^^(s). 

Several fits are performed, varying the minimum fitting 
scale in order to study the impact of non-linearities on the 
accuracy of the measurements as a function of scale. In all 
cases the maximum scale in the fitting process is limited to a 
given Smax, to avoid the complication of modelling the BAO 
signature. 

In the case of fitting the linear Kaiser/Hamilton model, 
there is only a single free parameter, /3. For the Dispersion 
model, we must also deal with a nuisance parameter in the 
form of the pairwise dispersion, ai 2 . The constraints on /3 
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Figure 9. Same as in Figure 8 but now fitting monopole and quadrupole moments of the two-point correlation function in redshift 
space. Continuous and dashed lines correspond respectively to the dispersion and pure linear models. Note how the linear model fails 
whenever galaxies are involved using either their auto-correlation or the group-galaxy cross-correlation, while it produces a meaningful 
result when applied to groups. 


that result from this model are marginalized by integrating 
the likelihood over a range of cri 2 , assuming a uniform prior 
within a maximum of <712 — 10/i“^Mpc. To some extent, 
there is a degeneracy between this parameter and (5: raising 
/3 battens the contours of whereas the Fingers of God 
oppose this tendency. But in practice this degeneracy is not 
severe, reflected in the fact that the errors on /3 are not much 
greater for the Dispersion model than for the linear model. 
The signihcance of treating Fingers of God is therefore more 
in helping to reduce bias in the best-htting value of (5. 


4.1 Fits to the full anisotropic correlation function 

We show here the results of htting the full tt). Fits are 

performed varying the minimum transverse scale and 

with a maximum scale hxed at — 50/i“^Mpc. 

We maximize the likelihood function dehned in equation 
(35), restricted to the diagonal elements of the covariance 
matrix. Following previous work Hawkins et al. (2003), 
Guzzo et al. (2008), Bianchi et al. (2012) we ht the quantity 

y{rp, tt) = ln[l + F{rp,Tr)] (42) 

rather than directly the values of ^^(rp,7r), in order to en¬ 
hance the weight of large, more linear scales. 

The results for the galaxy auto-correlation (blue lines), 
group-galaxy cross-correlation (green lines) and the group 
auto-correlation (red lines) are presented in Figure 8 and the 
related statistical errors are shown with the error bars. Gon- 
tinuous lines with hhed points correspond to hts performed 
using the Dispersion model while the dashed lines show the 
results when using the pure linear Kaiser/Hamilton model. 

When using the Dispersion model, the galaxy auto¬ 
correlation approach underestimates /3gai by ~ 10%, while 
the group-galaxy cross-correlation gives similar results when 
including scales below 10/i“^Mpc. The recovered values of 
/3gai from these two approaches are compatible with each 
other. On the other hand, the group auto-correlation results 
yield a much more accurate estimate of the group distortion 
parameter /3gr, underestimating it by 2-4%. As expected, 
given the different size of the data sets in these different 
cases, the statistical errors increase passing from the galaxy 


auto-correlation to the group-galaxy cross-correlation to the 
group auto-correlation. Results from the auto-correlations 
of galaxies and groups are in full agreement with the previ¬ 
ous work of Okumura & Jing (2011), de la Torre & Guzzo 
(2012) and Bianchi et al. (2012). In particular. Figure 5 in 
the last of these papers shows how the classical approach 
underestimates by ~ 10% for the case of haloes with mass 
~ 10^^Mq and that this systematic error diminishes for in¬ 
creasing halo masses, being close to unbiased for ~ 10^^Mq 
i.e. for group-sized haloes. 

On scales below 15 /i“^Mpc, the pure linear model heav¬ 
ily underestimates the distortion parameter with respect 
to the dispersion model when applied to the galaxy auto¬ 
correlation. Such discrepancy decreases using hrst the group- 
galaxy cross-correlation and disappears almost completely 
in the case of the group auto-correlation. This shows how 
the corrections to the linear model on intermediate quasi- 
linear scales, made through the dispersion model, become 
gradually less important using objects less affected by the 
non-linear effects on such scales. 

The results from the analysis in this section match 
our initial expectations. Specihcally, using objects tracing 
higher-mass haloes we are in general less sensitive to the 
details of non-linear corrections, as previously shown by 
Bianchi et al. (2012). This is particularly true for the auto¬ 
correlation of groups, where we see that the linear and Dis¬ 
persion models perform similarly even when including small 
scales. This is at variance with what we observe for galax¬ 
ies. The use of the cross-correlation, in this case, represents 
a reasonable compromise between having a larger statistics 
(thus smaller error bars), while limiting the systematic er¬ 
rors. In the next session, we shall see how a different way 
of htting the data can further ameliorate these results, in 
particular for the cross-correlation function. 


4.2 Fits to the standard multipole moments 

We perform here joint hts of the monopole (s) and the 
quadrupole (s) of the 2PGF in redshift space. We do 
not consider higher order moments, which are too noisy. To 
reduce scale dependence, rather than htting the multipoles 
directly, we consider the quantity (e.g. de la Torre et al. 
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Figure 10. Same as in Figure 8 but now fitting truncated monopole and quadrupole moments of the 2PCF in redshift space. 


2013) 

(43) 

which we fit between a varying minimum scale Smin and a 
maximum separation Smax = 80/i“^Mpc. Also in this case 
we limit the likelihood function in equation (35) to the diag¬ 
onal elements of the covariance matrix. As in the previous 
case, we fit using both the full Dispersion model and the 
simple linear model only, with results plotted as solid and 
dashed lines, respectively, in Figure 9. 

The measurements show a different behaviour, when 
compared to the fits to the full tt). Overall, the results 

are much more sensitive to the minimum fitting scale Smin- 
Also, the linear model gives highly biased results whenever 
galaxies are involved (either using their auto-correlation or 
the group-galaxy cross-correlation), producing a meaningful 
result only when applied to groups alone. 

The measurements from the galaxy auto-correlation 
(blue line) and the group-galaxy cross-correlation (green 
line) are now very similar, with an overall systematic er¬ 
ror on /3, which remains confined to values smaller than 
5%. On the other hand, the estimates from the group 
auto-correlation function have a highly scale-dependent be¬ 
haviour, with a significant positive bias when including 
scales below 10 /i“^Mpc. All approaches converge to a sys¬ 
tematic (negative) error of 3 — 5% when using only scales 
> 25 h~^Mpc. 

We note how in this case, compared to the analysis us¬ 
ing the full ^^(rp,7r), there is no clear indication that one 
method performs better than another. We interpret this as 
the consequence of the projection of the 2PCF onto the Leg¬ 
endre polynomials, which re-distributes over all scales s the 
non-linear effects originally mostly confined to small trans¬ 
verse scales Tp. As a consequence, limiting the fitting range 
to be above a given Smin does not eliminate such small-scale 
contribution. 

4.3 Fits to the truncated multipole moments 

As described in Section 2.2, we apply now the truncated mul¬ 
tipole moments (equation 28). These have been specifically 
defined as to eliminate the contribution of small transverse 
scales Tp, which, as we have just seen, affect all separations s. 
As in the previous section, we perform a joint fit of the quan¬ 
tity in equation (43) for the case of the truncated monopole 


£ = 0 and quadrupole £ — 2 using only diagonal errors and 
for the usual varying ranges in scale. 

The results, plotted in Figure 10, show a quite different 
situation from that of the two previous sections. The esti¬ 
mates using the galaxy auto-correlation function (blue lines) 
are not improved compared to using the full multipoles, but 
maintain a typical negative systematic error of ~ 5%. Con¬ 
versely, the estimates using the group auto-correlation (red 
lines) show a positive bias for any fitting range, converging 
to percent errors when including only scales larger than 15 
/i“^Mpc. Finally, in this case the fits to the cross-correlation 
are surprisingly stable, almost independently of the fitted 
range, with an expected systematic error of 3 % or less for 
Smin > 10/i“^Mpc. 

As evident from the dashed lines, using the linear model 
alone gives in general virtually identical results to the fit us¬ 
ing the dispersion model. The only marginal exception is the 
case of the galaxy auto-correlation. This overall behaviour 
indicates how the truncated multipoles are able to suppress 
the weight of small-scale non-linearities, thus making the 
role of the dispersion factor negligible. In the case of the 
galaxy auto-correlation, there is still a difference between 
the two approaches when including scales below 15 h~^Mpc 
in the fit, indicating the stronger effect of non-linearities in 
this case, compared to when groups are involved. Overall, 
we can conclude that the newly defined statistic of trun¬ 
cated multipole moments is helpful in reducing the impact 
of non-linearities on all scales in the measurements. 


4.4 Full covariance matrix analysis 

Measurements of the two-point correlation function in two 
different bins i and j are, in general, correlated with each 
other. Keeping in mind that the estimate and use of a proper 
covariance matrix is a non-trivial issue (see e.g. de la Torre 
& Guzzo 2012), we explore here the impact of including the 
full covariance matrix in the analysis of the standard and 
truncated multipole moments. Specifically, we estimate and 
use the joint covariance matrix of the truncated monopole 
{£ = 0) and quadrupole {£ = 2). This means measuring 
not only the covariance between the measurements of the 
monopole y®’^°^(s) and quadrupole y®’^^^(s) in two different 
bins separately but also the cross-covariance between the 
measurement of the monopole in a given bin i' and 

that of quadrupole (s) in bin j'. To do this we store the 
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Figure 11. Same as in Figure 9 but now fits are performed using the joint covariance matrix of standard monopole and quadrupole 
moments of the 2PCF in redshift space. 



Figure 12. Measurements (points with error-bars) and the related best fit models (lines) for the standard monopole (left panels) and 
quadrupole (right panels). Blue squares and dashed lines represent the auto-correlation of (HOD, L > IL*) galaxies, red triangles and 
dash-dotted lines represent the auto-correlation of (M > Mq dark matter haloes) groups while their cross-correlation is shown 

through green circles and continuous lines. Each row shows fit performed using the Dispersion model and full covariance matrix at 
different minimum fitting scales Smin- 


measurements of the monopole and the quadrupole 

concatenating them into a single array. The covari¬ 
ance matrix is then measured using the definition in equa¬ 
tion (34) from 27 sub-sample realizations. We recall that to 
avoid a singular covariance matrix and keep a good spatial 
resolution in our measurements we keep fixed the number of 
bins in the fitting range [smzn, 80] /i“^Mpc, independently of 

Smin • 

We repeat the analysis of standard and truncated mul- 
tipole moments of the 2PCF but now including the full data 
covariance matrix in the fitting procedure. The best-fitting 


models are presented in Figures 12 and 14, and the system¬ 
atic errors on /3 in Figures 11 and 13. In general, we find 
the behaviours with minimum fitting scale for the different 
types of correlation to be similar to the case where only di¬ 
agonal elements of the covariance matrix are used. The sta¬ 
tistical errors on /3 are however larger, as expected, and the 
detailed dependence on minimum fitting scale slightly nois¬ 
ier, consistently with the increased statistical errors. The 
only marginal difference in systematics compared to the 
diagonal covariance case, is in the galaxy auto-correlation 
where we note that both Kaiser and Dispersion models re- 
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Figure 13. Same as in Figure 10 but now fits are performed using the joint covariance matrix of truncated monopole and quadrupole 
moments of the PCF in redshift space. 




Figure 14. Same as in Figure 12 but for truncated multipole moments. 


cover lower values of 13 than previously, in particular for 
Smin ^ 20/i“^Mpc. This analysis confirm the results ob¬ 
tained previously and based on using diagonal covariance 
matrix. It confirms in particular the significantly improved 
performances of the truncated multipole moments of the 
galaxy-group cross-correlation with linear Kaiser model. 


4.5 Dependence of the results on the mass 
threshold of the group catalogue 

So far we have used as “groups” a catalogue of dark matter 
haloes with M > Mq. In this section we test how 

strong is the dependence of the main results obtained so 
far, on this mass threshold. To this end, we create the set 


of group catalogues listed in table 2. We limit our tests to 
using the truncated multipole moments of the 2PCF and 
obviously include the full covariance matrix of the data. The 
results for the group auto correlation are plotted in Figure 
15. These show a monotonic trend of the systematic error, 
almost independently of the minimum fitting scale Smin • The 
“sweet spot” for which the error is minimised appears to 
correspond to masses around h~^M q. This behaviour 

agrees with the previous result obtained by Bianchi et al. 
(2012) using halo catalogues from the BASICC simulation. 

A similar trend is seen in the group-galaxy cross corre¬ 
lation (Figure 16), but in this case there is a dependence on 
the minimum fitting scale and a general tendency to under¬ 
estimate the value of the distortion parameter for all group 
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Figure 15. Global dependence of the systematic errors on /3 obtained from the group 2PCF, when the mass threshold of the group 
catalogue is increased (Table 2). These tests use the truncated multipoles and include the full covariance matrix. As in previous plots, the 
dashed lines and empty markers in the left plot correspond to using the simple Kaiser model fit. The right panel is a zoom in, including 
only the Dispersion model curves, to avoid confusion. The red curve (M > Mq) corresponds to the group catalogue used so far 

in the analyses. 




Figure 16. Same as in Figure 15, but now for the cross-correlation of groups and galaxies as defined before. The values obtained from 
the galaxy 2PCF are also shown as blue squares and lines, for comparison. The colour coding for other curves is the same as in Figure 


15. 


catalogues. From this figure, we see that the group catalogue 
used for most of the tests in the paper (i.e. M > 10^^h~^M q 
- green line) is fairly representative of the general behaviour 
of the cross-correlation function when used to estimate (3. 


5 DISCUSSION AND CONCLUSIONS 

This work has explored two different ways of improving the 
accuracy of the measurement of the growth rate of cosmo¬ 
logical structure: by using the cross-correlation of individ¬ 
ual galaxies with groups of galaxies as well as by using a 
novel estimator of the two-point statistics in redshift space, 
the truncated multipole moments. The aim is to reduce the 
impact of non-linearities arising from small-scale random 
peculiar pairwise velocities, with respect to the usual ap¬ 
proach of using the multipole moments of the galaxy auto¬ 
correlation. We have used a set of simulated catalogues of 
galaxies to compare the accuracy with which the anisotropic 
2PCF tt), its multipole moments and its trun¬ 
cated multipole moments allow the recovery of the 

RSD parameters (3. In this comparison we compared both 
linear theory and the Dispersion model for RSD. 

We find that fitting the full anisotropic auto-correlation 
function of galaxies underestimates the distortion parame¬ 
ter by about 10%, confirming the results of Bianchi et al. 


(2012). The group-galaxy cross-correlation reduces this bias 
to a level of about 7-8%, and the group auto-correlation 
(for groups with mass larger than Mq) provides us 

with even less biased results, reaching an accuracy of about 
5%. As one may have expected, there is almost no difference 
between using either linear theory or the Dispersion model 
when fitting the group auto-correlation function: the Finger 
of God effect on group centroids is rather minor. 

The analysis of standard multipole moments gives no 
clear indications of the best choice of tracer for RSD. While 
the galaxy auto-correlation and the group-galaxy cross¬ 
correlation lead to similar results, the group auto-correlation 
produces highly scale-dependent measurements of the dis¬ 
tortion parameter. Such complications can mostly be ex¬ 
plained by the fact that strong non-linearities on small trans¬ 
verse scales Tp are integrated over all scales when projecting 
the 2PCF on Legendre polynomials, which are not captured 
by linear or dispersion RSD models. This does not happen 
when fitting the full ^(rp,7r) after excluding scales below a 
given rp. 

This problem is alleviated by the truncated multi¬ 
pole moments that we have introduced. These provide 
the most accurate estimates of the distortion parameter 
among the present tests, improving the accuracy over the 
use of ^^(rp,7r) and ^^’^^^(s). In that case, the galaxy 
auto-correlation underestimates the distortion parameter by 
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1^0g,[M>rnin / {h ^A7q)] 

N 

bgr 

ofid 

R>gr 

bi2 

12.00 

3,439,747 

1.041 ±0.001 

0.5220 ± 0.0005 

1.279 ± 0.004 

12.25 

1,997,835 

1.133 ±0.001 

0.4796 ± 0.0004 

1.175 ±0.004 

12.50 

1,140,387 

1.254 ±0.002 

0.4333 ± 0.0007 

1.061 ± 0.004 

12.75 

640, 598 

1.402 ± 0.002 

0.3876 ± 0.0006 

0.949 ± 0.003 

13.00 

350,518 

1.588 ± 0.003 

0.3422 ± 0.0006 

0.838 ± 0.003 

13.25 

184,238 

1.810 ±0.004 

0.3002 ± 0.0007 

0.735 ± 0.003 

13.50 

92,403 

2.094 ± 0.006 

0.2595 ± 0.0007 

0.636 ± 0.003 


Table 2. Main parameters for the catalogues of the dark matter haloes used as proxies for galaxy groups in Section 4.5. The first column 
contains the minimum mass limit for each catalogue, N is the number of objects, bgr is the linear bias, is the fiducial value for the 
distortion parameter and 612 is the relative bias between galaxies and groups as defined in Equation (23). 


about 5%, while the group auto-correlation becomes gradu¬ 
ally less biased when using larger minimum scales to reach a 
few percent accuracy on scales greater than 20 h~^Mpc. The 
group-galaxy cross-correlation, when similarly fitted with 
the Dispersion model, produces more stable and also more 
accurate measurements of reaching the percent level of ac¬ 
curacy when fitting scales greater than 15/i“^Mpc. A com¬ 
parison with the results from the linear model shows how 
the truncated multipole moments allow reducing the impact 
of small-scale non-linearities in RSD measurements, making 
possible the analysis with the simple linear model. In fact, 
for the galaxy auto-correlation, the limit at which the Dis¬ 
persion model breaks down is translated to 15 /i“^Mpc, com¬ 
pared with 20/i~^Mpc in the case of the anisotropic 2PCF. 

We studied the impact of bin-to-bin covariances in the 
fitting procedure and found no significant difference in terms 
of systematics, confirming our findings based on using only 
the diagonal covariance matrix. 

Finally, we have directly tested how our general results 
may depend on the mass threshold chosen to define the the 
group catalogue used for most of the tests 
The exercise using a set of further six group catalogues, 
with minimum mass thresholds ranging from 10 ^^h~^M q 
to h~^M q, shows that the “standard” group catalogue 
is fairly representative of the general trend observed in the 
systematic errors when using the cross correlation function. 
Additionally, it further confirms the dependence of the er¬ 
rors on the halo mass threshold, evidenced in Bianchi et al. 
( 2012 ). 

Although none of the methods studied here yields a zero 
bias, we find the results encouraging. Small systematics at 
the few percent level arise in RSD from other effects to do 
with the sky sampling (de la Torre et al. 2013), and these 
are already corrected for by analysis of mock data. The same 
approach could be taken in order to incorporate small sys¬ 
tematic errors in the theoretical RSD models being used, 
although further work will be required in order to demon¬ 
strate that these systematic offsets themselves are consistent 
independent of the true cosmology under study. 

The other element of this study that could benefit from 
extension concerns the group catalogue. The present work 
is somewhat idealised in that the group proxies are dark- 
matter haloes that are found directly in the simulation using 
more information than would be available in a real galaxy 
survey. The next step is therefore to repeat this analysis 
using a full simulation of the construction of an empirical 
group catalogue by linking the simulated galaxies in redshift 


space (following e.g. Robotham et al. 2011). At large group 
masses, uncertainties in group centroids should not be large, 
and so the features of the present work in terms of small 
Fingers of God should be reproduced. 
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